Incoherent Quantum Control 



Alexander Pechen* and Herschel Rabitz^ 



Department of Chemistry, Princeton University, 
Princeton, New Jersey 08544, USA 



Abstract 



Conventional approaches for controlling open quantum systems use coherent control 
which affects the system's evolution through the Hamiltonian part of the dynamics. 
Such control, although being extremely efficient for a large variety of problems, has 
limited capabilities, e.g., if the initial and desired target states have density matrices 
with different spectra or if a control field needs to be designed to optimally transfer 
different initial states to the same target state. Recent research works suggest extending 
coherent control by including active manipulation of the non-unitary (i.e., incoherent) 
part of the evolution. This paper summarizes recent results specifically for incoherent 
control by the environment (e.g., incoherent radiation or a gaseous medium) with a 
kinematic description of controllability and landscape analysis. 

1 Introduction 

The manipulation of atomic or molecular quantum dynamics commonly uses coherent quan- 
tum control, which may be extremely useful for a large variety of problems [1, 2, 3, 4, 5, 6, 
7, 8, 9, 10, 11, 12]. The dynamical evolution of a closed quantum system under the action 
of a collection of coherent controls u = {ui{t)} (e.g., Rabi frequencies of the applied laser 
field) is described by the equation 



Here pt is the system density matrix at time t (for an n-level quantum system, the set of all 
density matrices is = {p G Ain | P > 0, Tr p = 1}, where Ain = C"^" is the set of n x n 
complex matrices), Hq is the free system Hamiltonian describing evolution of the system in 
the absence of control fields and each Qi is an operator describing the coupling of the system 
to the control field ui{t). 




P\t=o = Po 
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Coherent control of a closed system induces a unitary transformation of the system 
density matrix pt = UtPoU} and may have some limitations. The first limitation is due 
to the fact that unitary transformations of an operator preserve its spectrum; thus the 
spectrum of pt is the same at any t and, for example, a mixed state po will always remain 
mixed [13]. A second limitation is that a control ?7,opt which is optimal for some initial state 
Po may be not optimal for another initial state po even if po and po have the same spectrum. 
This limitation originates from the reversibility of unitary evolution and is due to the fact 
that Utp^Ul ^ UtPoU} if Po ^ Po. To overcame these hmitations at least to some degree, 
control by measurements [14, 15. 16. 17, 18] or incoherent control [19, 20, 21, 22] may be 
used, and in this work incoherent control by the environment [19] (ICE) is discussed. Some 
general mathematical notions for the controlled quantum Markov dynamics are formulated 



The necessity to consider incoherent control relies also on the fact that coherent control of 
quantum systems (e.g., of chemical reactions) in the laboratory is often realized in a medium 
(solvent) which interacts with the controlled system and plays the role of the environment. 
Furthermore, then environment may be also affected to some degree by the coherent laser 
field, thus effectively realizing incoherent control of the system. Moreover, laser sources of 
coherent radiation at the present time have practical limitations, and some frequencies are 
very expensive to generate compared to the respective sources of incoherent control (e.g., 
incoherent radiation as considered in Sec. 2.1 of this work). Thus the latter incoherent 
control can be used in some cases to reduce the total cost of quantum control. 

This paper summarizes recent results specifically for incoherent control by the environment [19] 
(ICE). A general theoretical formulation for incoherent control is provided in Sec. 2, followed 
by the examples of control by incoherent radiation (Sec. 2.1) and control through collisions 
with particles of a medium (e.g., solvent, gas, etc.. Sec. 2.2). Relevant known results about 
controllability and the structure of control landscapes for open quantum systems in the 
kinematic picture are briefly outlined in Sec. 3. 

2 Incoherent control by the environment 

The dynamical evolution of an open quantum system under the action of coherent controls 
in the Markovian regime is described by a master equation 



The interaction with the environment modifies the Hamiltonian part of the dynamics by 
adding an effective Hamiltonian term ifeff to the free Hamiltonian Hq. Another important 
effect of the environment is the appearance of the term C which describes non-unitary aspects 
of the evolution and is responsible for decoherence. This term in the Markovian regime has 
the general Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) form [24, 25] 



in Ref. [23]. 
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where Lj are some operators acting in the system Hilbert space. The exphcit form of the 
GKSL term depends on the particular type of the environment, on the details of the mi- 
croscopic interaction between the system and the environment, and on the state of the 
environment. 

The coherent portion of the control in (2) addresses only the Hamiltonian part of the evo- 
lution while the GKSL part L remains fixed (for the analysis of controllability properties for 
Markovian master equations under coherent controls see for example, Ref. [26]). However, 
the generator L can also be controlled to some degree. For a fixed system-environmental 
interaction, the generator H depends on the state of the environment, which can be either a 
thermal state at some temperature (including the zero temperature vacuum state) or an ar- 
bitrary non-equilibrium state. Such a state is characterized by a (possibly, time dependent) 
distribution of particles of the environment over their degrees of freedom, which are typi- 
cally the momentum k e and the internal energy levels parameterized by some discrete 
index a ^ A (e.g., for photons a = 1, 2 denotes polarization, for a gas of A^-lcvel particles 
a = 1,...,N denotes the internal energy levels). Denoting the density at time t of the 
environmental particles with momentum k and occupying an internal level a by nk,a(t), and 
the corresponding GKSL generator as £ = £[nk,a(t)], the equation (2) becomes 



Here both ui{t) and nk,a(t) are used as the controls, and for ny_^a{t) the optimization is done 
over k, a in a time dependent fashion to obtain a desired outcome. 

The solution of (3) with the initial condition p\t=o = Po for each choice of controls {ui{t)} 
and nk,a(^) can be represented by a family Pt{{ui),nk^a}, ^ > of completely positive (CP), 
trace preserving maps (see Sec. 3 for the explicit definitions) as 



In general, for time dependent controls this family forms not a semigroup but a self-consistent 
two-parameter family of CP, trace preserving maps Pt,T{iui),nk^a}, t > t > 0, where each 
Pt,T{{ui),n]i,a} represents the evolution from r to t. 

The target functional, also called the performance index, describes a property of the 
controlled system which should be minimized during the control and commonly consists of 
the two terms: 



The term Ji[(ii/), nk,a], called the objective functional, represents the physical system's prop- 
erty which we want to minimize. The term J2[{ui),nk^a\, called the cost functional, represents 
the penalty for the control fields. 

The first general class of objective functionals appears in the problem of minimizing 
the expectation value of some observable associated to the system at a target time T > 0. 
The system is assumed at the initial time t = to be in the state po- Any observable 
characterizing the system (e.g., its energy, population of some level, etc.) is represented 
by some self-adjoint operator O acting in the system Hilbert space, and the corresponding 
objective functional has the form 




(3) 



Pt = Pt{{ui),nk,a}Po 



(4) 



J[{ui),nk,a] = -^l[(lii)>"'k,a] + -^2[(ti/),?T'k,a] 



-^i[(iii)> ^k,a] = Tr [pt{ui, nk,a)0] = Tr [{PT{{ui),nk^a}po)0] 



(5) 
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Here PT{ui,nk^a) = -Pr{('Ui)> ^k,a}po is the final density matrix of tlie system evolving from 
the initial state po under the action of ui and nk,a. The physical meaning of this objective 
functional is that it represents the average measured value of the observable O at the final 
time T when the system evolves from the initial state po under the action of the controls 

The second general class of objective functionals appears in the problem of optimal state- 
to-state transfer. Suppose that initially the system is in a state po and the control goal is 
to steer the system at some target time T into some desired target state ptarget- In this case 
one seeks controls (ui) and nk,a: which minimize the distance between the states Pt{ui, 11^,0} 
and ptarget- The corresponding objective functional has the form 

-^l[K)>"'k,a] = \\PT(ui,nk^a) " PtargetH = \\PT{(ui),nk^a}Po " PtargetH (6) 

where || ■ || is a suitable matrix norm. Usually the Hilbert-Schmidt norm \\A\\ = Vr^AfA 
can be used. 

The third important class of objective functionals appears in the problem of producing 
a desired target CP, trace preserving map Ptarget- In this case the objective functional has 
the form 

-^l[(«0>^k,a] = ||Pr{(Mz),nk,a} - Ptargetll (7) 

where || • || is a suitable norm in the space of all CP, trace preserving maps. In particular, in 
conventional models of quantum computation the target transformation Ptarget is a unitary 
gate (e.g., a phase or Hadamard C/h gate, and for these examples Ptarget = U^f, or Ptarget — 
Um, respectively) [27, 28]. More general non-unitary target transformations can arise [e.g., 
in quantum computing with mixed states [29] or for generating controls robust to variations 
of the initial system's state [30] (sec also Sec. 3.1)]. 

The cost functional J2 can be chosen to have the form 

-^2[(ii/),rik,i] = / dtai{t)\ui{t)f + max / (ik/3j(k)nk,j(t) 
I Jo i J 

Here each function ai{t) > [resp., /?i(k) > 0] is a weight describing the cost for the control 
ui at time t (resp., for the density of particles of the environment with momentum k and 
occupying the internal level i) . The first term minimizes the energy of the optimal coherent 
control. The second term minimizes the total density of the environment. 

The control functions belong to some sets of admissible controls (ui) e £ and nk,a G V. 
The following three important problems arise. 

Optimal controls. Find, for a given initial state po and a target time T, some (or all) 
controls (ui) e £ and nk,a £ ^ which minimize the performance index. 

Reachable sets. Find, for a given final time T > and an initial state po, the set of all 
states reachable from po up to the time T, i.e., the set 

'^r(po) = {Pt{{ui),nk,a}po \t<T, (ui) e £, nk,a e T>} 

Landscape analysis. Find, for a given T > 0, an initial state po and a self-adjoint op- 
erator O, all extrema (global and local, and saddles, if any) of the objective functional 
Ji[ui(t) , n^^ait)] = Tr [{Pr [(«;), n]po}0] defined by (5) [and similarly for the objective func- 
tionals defined by (6) and (7)] or of the corresponding performance index. 
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2.1 Incoherent control by radiation 

Non-equilibrium radiation is characterized by its distribution in photon momenta and polar- 
ization. For control with distribution of incoherent radiation the magnitude of the photon 
momentum |k| can be exploited along with the polarization and the propagation direction in 
cases where polarization dependence or spatial anisotropy is important (e.g., for controlling 
a system consisting of oriented molecules bound to a surface). 

A thermal equilibrium distribution for photons at temperature T is characterized by 
Planck's distribution 

1 

nk = - 



( ) - 1 



where c is the speed of light, h and k-Q are the Planck and the Boltzmann constants which we 
set to one below. Non-equilibrium incoherent radiation may have a distribution given as an 
arbitrary non-negative function nk,a(i). Some practical means to produce non-equilibrium 
distributions in the laboratory may be based either on filtering thermal radiation or on the 
use of independent monochromatic sources. 

The master equation for an atom or a molecule interacting with a coherent electromag- 
netic field Ec{t) and with incoherent radiation with a distribution nk(t) in the Markovian 
regime has the form: 

^ = -i[H^ + - f^E,{t),pt] + Cji^M)]pt (8) 

The coherent part of the dynamics is generated by the free system's Hamiltonian Hq = 
'^n^nPn with eigenvalues £„, forming the spectrum spec(ifo), and the corresponding pro- 
jectors Pn, the effective Hamiltonian H^s resulting from the interaction between the system 
and the incoherent radiation, dipole moment fi, and electromagnetic field Ec{t). 

The GKSL generator C = C^ad induced by the incoherent radiation with distribution 
function nk(t) has the form (e.g., see Ref. [31]) 

'^^Rad[^k(i)]p = + 7la;(^)](2/ia;P/^t - /J-ll^u^P - P/J-ll^uj) (9) 

Here the sum is taken over the set of all system transition frequencies = {£„ — £^ | £„, e 
specifo}, A*a; = Pm/J'Pn, ^nd the coefficients 



l^(t) = 7rj dkS{\k\-oj)\gi,\^[nk{t) + (1 ± l)/2] 

determine the transition rates between energy levels with transition frequency lv. The tran- 
sition rates depend on the photon density nk(t). The form- factor (^k determines the coupling 
of the system to the k-th mode of the radiation. Equation (8) together with the explicit 
structure (9) of the GKSL generator provides the theoretical formulation for analysis of 
control by incoherent radiation. 

The numerical simulations illustrating the capabilities of learning control by incoherent 
radiation to prepare prespecified mixed states from a pure state is available [19] along with a 
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theoretical analysis of the set of stationary states for the generator ^Rad for some models [21] . 
Incoherent control by radiation can extend the capabilities of coherent control by exciting 
transitions between the system's energy levels for which laser sources are either unavailable 
at the present time or very expensive compared with the corresponding sources of incoherent 
radiation. Ref. [22] provides a simple experimental realization of the combined coherent (by 
a laser) and incoherent (by incoherent radiation emitted by a gas-discharge lamp) control of 
certain excitations in Kr atoms. 

2.2 Incoherent control by a gaseous medium 

This section considers incoherent control of quantum systems through collisions with particles 
of a surrounding medium (e.g., a gas or solvent of electrons, atoms or molecules, etc.). 
This case also includes coherent control of chemical reactions in solvents if the coherent 
field addresses not only the controlled system but the solvent as well. The particles of the 
medium in this treatment serve as the control and the explicit characteristic of the medium 
exploited to minimize the performance index is in general a time dependent distribution 
of the medium particles over their momenta k and internal energy levels a & A. This 
distribution is formally described by a non- negative function n : x ^4 x R — > R+, whose 
value nk,a(^) (where k e R^, a 6 A, and t e R+) has the physical meaning of the density at 
time t of particles of the surrounding medium with momentum k and in internal energy level 
a. In this scheme one prepares a suitable, in general non-equilibrium, distribution of the 
particles in the medium such that the medium drives the system evolution through collisions 
in a desired way. 

It may be difficult to practically create a desired non-equilibrium distribution of medium 
particles over their momenta. In contrast, a non-equilibrium distribution in the internal 
energy levels can be relatively easily created, e.g., by lasers capable of exciting the internal 
levels of the medium particles or through an electric discharge. Then the medium particles 
can affect the controlled system through collisions and this influence will typically depend 
on their distribution. A well known example of such control is the preparation of population 
inversion in a He-Ne gas-discharge laser. In this system an electric discharge passes through 
the He-Ne gas and brings the He atoms into a non-equilibrium state of their internal degrees 
of freedom. Then He-Ne collisions transfer the energy of the non-equihbrium state of the 
He atoms into the high energy levels of the Ne atoms. This process creates a population 
inversion in the Ne atoms and subsequent lasing. A steady electric discharge can be used 
to keep the gas of helium atoms in a non-equilibrium state to produce a CW He-Ne laser. 
This process can serve as an example of incoherent control through collisions by considering 
the gas of He atoms as the control environment (medium) and the Ne atoms as the system 
which we want to steer to a desired (excited) state. 

Quantum systems controlled through collisions with gas or medium particles in certain 
regimes can be described by master equations with GKSL generators whose explicit structure 
is different from the generator >CRad describing control by incoherent radiation. If the medium 
is sufficiently dilute, such that the probability of simultaneous interaction of the control 
system with two or more particles of the medium is negligible, then the reduced dynamics of 
the system will be Markovian [32, 33] and will be determined by two body scattering events 
between the system and one particle of the medium. Below we provide a formulation for 
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control of quantum systems by a dilute medium, although the assumption of diluteness is 
not a restriction for ICE, and dense mediums might be used for control as well. 

The master equation for a system interacting with coherent fields ui{t) and with a di- 
lute medium of particles with mass m has the form (3) with the generator £[nk ,^(t)] = 
-^Medium [^k,a(^)] Specified by the distribution function of the medium ny_^a{t) and by the T- 
operator (transition matrix) for the scattering of the system and a medium particle. Below 
we assume that the particles of the medium are characterized only by their momenta and 
do not have internal degrees of freedom; otherwise, the state of one particle of the medium 
should have the form |k, a), where a specifies the internal degrees of freedom. A transition 
matrix element is T„.„/(k, k') = (n, k|T|n', k'), where |n,k) = |T7.)|k) denotes the product 
state of the system discrete eigenstate \n) (an eigenstate of the system's free Hamiltonian 
Hq with eigenvalue £„) and a translational state of the system and a medium particle with 
relative momentum k. If the system is fixed in space (we consider this case below correspond- 
ing to the system particle being much more massive than the particles of the surrounding 
medium) then |k) is a translation state of a medium particle. The general case of relative 
system medium particle motion can be considered as well using suitable master equations. 
We will use the notation T^(k, k') := Tm,n(k, k')\m){n\. The density of particles 

m,n: em-en=i^ 

of the medium at momentum k is denoted as rikit), and the set of all transition frequencies 
uj of the system among the energy levels of is denoted as ^2. In this notation the GKSL 
generator is 



UJ 



/■ r /Ik'P IkP 

^Medmm[nk{t)]p = Stt 2^ ^ dknk{t) J dk'5 - ^ + 

X "7L(k',k)prJ;(k',k) - ^{ri;(k',k)7L(k',k),p}] (lo) 

where {■, ■} denotes the anti-commutator. If the medium is at equilibrium with inverse tem- 
perature (3, then the density has the stationary Boltzmann form nk(t) = = C{P, n) exp(— /3|kp/2m). 
Here the normalization constant C(/3, n) is determined by the condition J dkrik = n, where n 
is the total density of the medium. The structure of Eq. (10) has been discussed previously for 
equilibrium media [32, 33] and for non-equilibrium stationary media [34]. Non-equilibrium 
media may be characterized by generally time dependent distributions. Equation (3) with 
>C[nk(t)] = >CMedium[^k(^)] providcs the general formulation for theoretical analysis of control 
by a coherent field ui{t) and by a non-equilibrium medium with density nk(t). 

As a simple illustration of such incoherent control. Fig. 1 reproduces the numerical results 
from Ref. [19] for optimally controlled transfer of a pure initial state of a four-level system 
into three different mixed target states [i.e., the objective function (6) is chosen]. The control 
is modelled by collisions with a medium prepared in a static non-equihbrium distribution n|k| 
whose form is optimized by learning control using a genetic algorithm (GA) [35] based on the 
mutation and crossover operations. Since the initial and target states have different spectra, 
they can not be connected by a unitary evolution induced by coherent control. However, 
Fig. 1 shows that ICE through collisions can work perfectly for such situations. 
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Figure 1: (Prom Ref. [19]. Copyright (2006) by the American Physical Society.) Results of 
ICE simulations with a surrounding non-equilibrium medium as the control for target states 
(a) ptargct = diag(0.3; 0.3; 0.2; 0.2), (b) ptargct = diag(0.3; 0.2; 0.3; 0.2), and (c) ptargct = 
diag(0.4; 0.1; 0.4; 0.1). Each case shows: the objective function vs GA generation, the 
optimal distribution vs momentum, and the evolution of the diagonal elements of the density 
matrix for the optimal distribution. In the plots for the objective function the upper curve 
is the average value for the objective function and the lower one is the best value in each 
generation. 



3 Kinematic description of incoherent control 

Physically admissible evolutions of an n-level quantum system can be represented by CP, 
trace preserving maps {Kraus maps) [36]. A map $ : A^„ Ain is positive if for any p G A4n 
such that p > 0: $(p) > 0. A linear map $ : Ain is CP if for any / G N the map 

$ ® I; : Mn ® M-i — M-n ® -Ml is positive (here I; denotes the identity map in Mi). A CP 
map $ is called trace preserving if for any p e Mn '■ Tr $(p) = Tr p. The conditions of trace 
preservation and positivity for physically admissible evolutions are necessary to guarantee 
that $ maps states into states. The condition of complete positivity has the following 
meaning. Consider the elements of Mi as operators of some /-level ancilla system which 
does not evolve, i.e., its evolution is represented by the identity mapping I/. Suppose that 
the n-level system does not interact with the ancilla. Then the combined evolution of the 
total system will be represented by the map $ (8) 1; and the condition of complete positivity 
requires that for any I this map should transform all states of the combined system into 
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states, i.e. to be positive. 

Any CP, trace preserving map $ can be expressed using the Kraus operator-sum repre- 
sentation as 

A 

1=1 

where Ki e (C"xri. Kraus operators subject to the constraint Yl!i=i ^l^i — to 

guarantee trace preservation. This constraint determines a complex Stiefel manifold T4(C'^") 
whose points are n x (An) matrices V = {Ki, K2; . . . ; Kx) (i.e., each y is a column matrix 
of Ki, . . . Kx) satisfying the orthogonality condition V'^V = I„. 

The explicit evolution -Pt{(M/), nk,a} in (4) is unlikely to be known for realistic systems. 
However, since this evolution is always a CP, trace preserving map, it can be represented in 
the Kraus form 

A 

Pt{{ui),nk^a}p = ^ Ki{t, {ui),nk^a)pKl{t, (i*;), nk,a), 
1=1 

A 

where ^Kj{t, {ui),nk^a)Ki{t, (Mz),nk,a}) = In 

i=l 

Assume that any Kraus map can be generated in this way using the available coherent 
and incoherent controls {ui{t)} and nk,a(^)- Then effectively the Kraus operators can be 
considered as the controls [instead of {ui{t)} and nk,a(^)] which can be optimized to drive 
the evolution of the system in a desired direction. This picture is called the kinematic picture 
in contrast with the dynamical picture of Sec. 2. In the next two subsections we briefly outline 
the controllability and landscape properties in the kinematic picture. 

3.1 Controllability 

Any classical or quantum system at a given time is completely characterized by its state. The 
related notion of state controllability refers to the ability to steer the system from any initial 
state to any final state, either at a given time or asymptotically as time goes to infinity, and 

the important problem in control analysis is to establish the degree of state controllability 
for a given control system. Assuming for some finite-level system that the set of admissible 
dynamical controls generates arbitrary Kraus type evolution, the following theorem implies 
then that the system is completely state controllable. 

Theorem 1 For any state pf G Vn of an n-level quantum system there exists a Kraus map 
$pj such that $pj (p) = Pf for all states p G I^„. 

Proof. Consider the spectral decomposition of the final state pf = XlILi -Pil*^*) where pi 
is the probability to find the system in the state (pi > and XlILi?** ~ -'-)• Choose an 
arbitrary orthonormal basis in the system Hilbert space and define the operators 
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The operators Kij satisfy the normahzation condition '^^j^i KjjKij = I„ and thus determine 
the Kraus map (p) = S"j=i KijpKjj. The map acts on any state p G Vn as 

n n 

*Pf(p) = ^PM{Xj\p\Xj){<l>i\ = (Trp) =pf 

i,j=l 1=1 

and thus satisfies the condition of the Theorem. □ 

The potential importance of this result is that it shows that there may exist a single 
incoherent evolution which is capable for transferring all initial states into a given target 
state, and moreover, the target state can be an arbitrary pure or a mixed state[30]. Thus 
this theorem shows that non-unitary evolution can break the two general limitations for 
coherent unitary control described in the second paragraph in the Introduction. 

3.2 Control landscape structure 

In the kinematic description, under the assumption that any Kraus map can be generated, the 
objective functional becomes a function on the Stiefel manifold T4(C'^"). In practice, various 

gradient methods may be used to minimize such an objective function. If the objective 
function has a local minimum then gradient based optimization methods can be trapped in 
this minimum and will not provide a true solution to the problem. For such an objective 
function, if the algorithm stops in some minimum one can not be sure that this minimum 
is global and therefore this solution may be not satisfactory. This difficulty does not exist 
if a priori information about absence of local minima for the objective function is available 
as provided by the following theorem for a general class of objective functions of the form 
Ji[Ki, . . . , Kx] = Tr [(Xl^=i ^iP^i)C^] i'^ the kinematic picture. 

Theorem 2 For any n eN, p E Vn, and for any Hermitian O G A4n the objective function 
Ji[Ki, . . . , Kx] = TT[{^'^^^KjpKj)0] on the Stiefel manifold Vn{C^"') does not have local 
minima or maxima; it has global minimum manifold, global maximum manifold, and possibly 
saddles whose number and the explicit structure depend on the degeneracies of p and O. 

The case n — 2 has been considered in detail in Ref. [37], where the global minimum, 
maximum, and saddle manifolds arc explicitly described for each type of initial state p. In 
particular, it is found that the objective function Ji for a non-degenerate target operator O 
and for a pure p (i.e., such that = p) does not have saddle manifolds; for the completely 
mixed initial state p = ^I, Ji has one saddle manifold with the value of the objective 
function Jgaddie = 1/2; and for any partially mixed initial state Ji has two saddle manifolds 
corresponding to the values of the objective function Jg^ddie ~ where w = Tr [per] 

and cr = (cTx, (TyjCTz) is the vector of Pauli matrices (the vector w is in the unit ball, ||w|| < 1 
and this vector characterizes the initial state as p = |[I2 -h (w, cr)]). The case of arbitrary n 
is considered in Ref. [38]. 

4 Conclusions 

This paper outlines recent results for incoherent control of quantum systems through their 
interaction with an environment. A general formulation for incoherent control through GKSL 
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dynamics is given, followed by examples of incoherent radiation and a gaseous medium 
serving as the incoherent control environments. The relevant known results on controllability 
of open quantum systems subject to arbitrary Kraus type dynamics, as well as properties of 
the corresponding control landscapes, are also discussed. 
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